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The semiclassical approximation for the Hamiltonian of Dirac particles interacting with an ar- 
bitrary gravitational field is investigated. The time dependence of the metrics leads to new con- 
tributions to the in-band energy operator in comparison to previous works on the static case. In 
particular we find a new coupling term between the linear momentum and the spin, as well as 
couplings which contribute to the breaking of the particle - antiparticle symmetry. 
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INTRODUCTION 



o 

t In this paper we consider the theory of Dirac fermions in an arbitrary curved space-time in the Hamiltonian 
• formulation. To reveal the physical content of the theory it is necessary to perform the diagonalization of the 
^-j. \ Hamiltonian uncoupling the positive and the negative energy states. For a fermion interacting with an electromagnetic 
— ■ field the Foldy-Wouthuysen (FW) transformation based on an approximate scheme valid in the non relativistic limit 
is often used [l| . This same method was also applied in all the previous studies of Dirac fermions in a gravitational 
field [2]. Here instead, we will consider the fully relativistic regime but in a semiclassical approximation for which 
the de Broglie wave length of the fermion must be much smaller that the characteristic size of the inhomogeneities of 
the external field. A recent semiclassical FW-like transformation used for Dirac particle in a strong electromagnetic 
field could be adapted to the gravitational problem 3]. But instead we will use another method developed by the 
authors which essentially differs from the FW. This method allows us to find the diagonal representation of any kind 
of matrix valued quantum Hamiltonian as a series expansion in the Planck constant. Here we will directly apply the 
general formula obtained at the semiclassical (first order in the Planck constant) limit to the case of Dirac fermions 
in an arbitrary curved spacetime. This is an extension of previous papers where massless and massive particles in a 
static gravitational fields were treated. The extension to time dependent metrics turns out to be non-trivial and leads 
to new coupling terms in the in-band energy operator which break the particle- antiparticle symmetry. 



ELECTRON IN A GRAVITATIONAL FIELD 

A one half spinning particle of mass m coupled to an arbitrary gravitational field is described by 4— spinor field ip 
satisfying the covariant Dirac equation 

(ih-f a D a - m)V = (1) 

sh ; 

where we use c = 1, but keep explicit the Planck constant h and a = 0, 1, 2, 3. 

The covariant spinor derivative is denned as D a — h l a Di with Di = di — \ [7°, 7^] r? . The matrices 7" are the 
usual Dirac matrices, h l a are the the orthonormal vierbein and T"^ the spin connection components. 



Rewriting Eq. (JXJ) under the Schrodinger form 



U&. = Hip 



Of 

we obtain the following Hamiltonian 

H = g m h%^ (rP a + m) + ~e B „ 7 rf £ 7 + ijr£% (2) 

where we introduced the notation for the pseudo-momentum P a =K b a (Pi + |e ei a 7 r^ /3 S 7 ) with the spin matrix X 7 , 
satisfying the relation e 7 "^S 7 = |(7 a 7 /3 — 7 /3 7 a )- We use the conventions of Bjorken and Drell Q for the Dirac 
matrices 7", and a? . 



2 



Surprisingly, Eq. ([2]) turns out to be non-hermitian for a time dependent metric. We then follow the approach of 
Leclerc [B| who showed that one must add the term |<9 t In (— gg 00 ) to make it Hermitian. It will be shown later on, 
that the presence of this term is also necessary for the diagonalization procedure to work. 

Therefore the Hamiltonian considered in the following will be 

H = goohh (l a P a +m) + \e ePl Tf^ + i^r° V + \d t In (-gg™) (3) 

The goal of this paper is therefore to diagonalise Eq. ([3]) to first order in h. But before embarking into this, we need 
first to discuss the definitions and properties of the scalar product in a curved space-time. 



Scalar product. 

As said before the Hamiltonian H is Hermitian. However, the notion of hermiticity is here defined with respect to 
a scalar product in curved space [5j , namely : 



(Mt) I ^2 (t)) m = J tPtV^9h° ^^ 2 d 3 x = I i;+U^ 2 d 3 x (4) 

where we introduced the notation U = ^/— 5/1^7^7 and tpi, ^2 two spinors. Therefore an operator O is hermitian 
for (Dy/^ defined in Eq. (g]), if 

itfV^W (Oik) = f (0^) + \/=^ 7 V (tfa) (5) 



It means that matricially 

+ U = UO (6) 

where "+" denotes from now, the usual Hermitic conjugation (transposition and complex conjugation), that is, the 
hermitic conjugate with respect to the scalar product in flat space denoted (|) and defined by 

tyi(t)liM*)>= / i>ti>2d z x (7) 



Unfortunatly the definition Eq. ([5]) turns out to be untractable for pratical computations. Actually, for the sake of the 
diagonalization procedure, we aim at working with matrices which are Hermitian with respect to the usual transpose 
and complex conjugate operation Eq. ([?]) ■ so that the diagonalization can be performed through a unitary matrix in 
the usual sense. To do so, notice that if O is Hermitian for Eq. (|3]), then Eq. © implies that U?OU~z is Hermitian 
for Eq. (0- Thus, starting with the Hamiltonian H defined in Eq. (|3|), U?HU~? is Hermitian in the usual sense and 
can be diagonalized through a standard unitary matrix (that is unitary for |7|). 

The Hamiltonian of interest for us will thus be U^HU^?. It's hermiticity in the usual sense allows us to write: 



H + 



U*,H 



u-i + \u- 



(8) 



f 

2 

The non unitarity of the transformation U is not problematic here, since it is precisely used to change the metric from 
curved to flat scalar product, and moreover H and U^HU^? have the same spectrum. In the case of a static metric 
(time independent) and satisfying = / (R) <5°, for a certain position dependent function / (R), the transformation 

t/2 reduces to the multiplication by a function of R. Then, using Eq. ([3]) for H, Eq. ([SJ simplifies easily to: 

1- 1 1 1 - 1 1 / 1 ~ i\+ I/~ 
U2HU-2 = -U*HU~* + -(UzHU-2) =-[H + H^ 

It is this form that was considered in If hi addition the metrics is diagonal, one recovers the transformation 

studied in Q. 

Independently of the practical advantages of the flat scalar product, there is an other and deeper reason to transform 
the Hamiltonian to the flat space. Actually, when diagonalizing the Hamiltonian with respect to Eq. ([7|) the diagonal 
subspaces of up and down spinors will appear to be obviously orthogonal. This is of course not the case for the scalar 
product Eq. (jj) which mixes both subspaces. As a consequence diagonalizing with respect to the curved space scalar 
product does not lead to a clear separation between particles and antiparticles. 
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Unitarity 

We will end up this section by stressing the fact that for a non-static metric the Hamiltonian and the time evolution 
operator do not coincide in the flat representation, and in addition, the time evolution operator cannot be made 
Hermitian. To make this point clearer, consider the Schrodinger equation 

d 

ih—^ = H^ 
dt 

Applying the transformation yields, 

JU' = (ui (<) ru-\ {t) - mu L * (t) lu-$ (t)\ & 

with *' = l/ity. As a consequence the evolution for 'J/', is given by the operator 

H e = C/5 (t) HU-i (t) - ihU? (t) (t) (9) 

One would like H e to be hermitian for the scalar product Eq. ([7]). However, due to the non unitarity of Lfi, the 
contribution ihU^ (t) -§[U~^ (t) is not. The reason for this non unitarity tracks back to the dependence in t of the 
scalar product Eq. (jj), so that the norm of a wave function is not preserved in time. Actually starting with an initial 
condition *' (to), the solution for the Schroedinger equation is 



ty'fa) = Texp Qf 1 H e (t) dtj *'(t ) = Ui (ti)Texp ^ ' 



ti 

H(t)dt]U-*(to)*' (t ) 



where T is the time ordered exponential. Assuming the norm *' (to) to be equal to 1, *' (ti) is easily seen to have a 
norm (respectively to Eq. (Q}) different from one. That can be checked easily on an infinitesimal timeslice, t\ = to+At. 
Indeed 



(*' (ti)| |*' (tx)) = (U* (ir) exp (-iH (to) At) * (t ) | (h) exp (iH (t Q ) At) * (t ) / 
= (exp (-iH (to) At) # (t ) | U (h) exp (iH (t ) At) * (t )) 

= (exp (-iH (to) At) * (t ) | exp (iH (t ) At) * (t )) u(tl) (10) 

and this is different from 1 since, 

(exp (-iH (t ) At) * (t ) | exp (iH (to) At) * (t )} u(tl) 

= (exp (-iH (t ) At) * (t ) | U (t ) exp (iH (t ) At) * (to)) + (t ) \ ^-U (t ) At* (t )J (11) 

= (exp (-iH (t ) At) * (t ) | exp (iH (t ) At) * (t )) u(to) + (* (to) \ (t ) At* (t )\ (12) 

the first term is equal to 1, actually * (to) is of norm 1 for (\)jj/ to -\ and exp (iH (to) At) is unitary for this scalar product. 
As a consequence, (*' (ti)| I*' (ti)) differs from one and H e is non unitary. The reason is clear from Eq. (fl0|) : a 
vector of norm 1 for (|)j/( to ) is transported to a vector, that has no more norm 1 for (|)y( tl v During the evolution, the 

matrix U defining the scalar product has changed too, and the non hermitian connexion term —ihlji (t) -j^U^^ (t) 
tracks the change of metric between t and t + At. 

Therefore, the time evolution of the state is non-unitary. In the rest of the paper we focus on the diagonalization of 
the energy operator U' (t) HU~' (t) this one being Hermitian, although the diagonalization of —ihU? (t) -§[U~2 (t) 
is provided for the sake of completness in appendix B. 



Transformation to the flat space 

We thus now focus on the Hamiltonian Eq. (j3|), and compute explicitly UiHU~i which as shown is hermitian in 
the usual sense. The transformation U is given by 

U = V=g (h° + h°a?) = y/=gh° (l + o^/fcg) 
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One can then deduce 



Ui =/(l+V) 



and thus 



where v? — upu^ with 



up 



and 



-9K 



„ 1+ v w ^ 



1/2 



V 



The greek (lorentzian) indices are assumed now to run only from 1 to 3. As a consequence the Hamiltonian H 
U^HU~^ given by Eq. §8§ reads 



which can be rewritten as 



H 



1 , , , 1 (l + Una 13 ) r fl1 1 r _l fli (1 + una 13 ) 

H = -(H + H+)- ^^-^ [H,u^] + - lH+,u^} ^-^> 



1 



/(1-u 2 ) 
1 



i-(l-uX) (Vp^ + ).v r — 



(1 - ^c/) 
(H-u^c/) 



/(1-tt 2 ) 

where P and R are the canonical momentum and position operator satisfying [i? l ,Pj] = ihSj. The Hamiltonian can 
also be written as 



H = l 



2 \(l-u 2 ) 



1 + 1 

H + H + 



1 1 



(l-u 2 )J 2(l-u 2 ) 



[H,upa p ] +i [ff+u^] - 1 



(1-u 2 ) 



H + H A 



l-u 2 



UpOi h 



l 2 



w/3a /3 ,(V P i?).V R 



(1 - ti 2 ) 
1 



/(l-u 2 ) 



(13) 



At this level, the computation of the last expression turns out to be quite technical is fully developped in Appendix 
A. The result is given by the following expression 



H 



1 



=-0/9 



r 



e/3 



r\i> 



1 f-n \ (goo (h° xh')-ux H 4 ) 
Jm + -g 1 ^ + P i9 * + h(r°+ T e j .£ + -— ^ ^. (Viu) J 



(14) 



where from now on, all indices i, f3,p... are only spatial and run from 1 to 3, but roman indices are raised and lowered 
by the meric gij and greek indices by the lorentzian metrics r) a p . The several notations introduced above are given by 
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the following expressions 

2g 00 (h° s hi - hill's) u s 

Tjl _ Tjl J V / 

ff = Tf - (H- 1 ) 1 ? gooh^ 



ff = (H- 1 )" H>t* p 



L j - ^[fePi 1 o + 9oo9 e QM l i + £ u 1 n p l 
^ 2 ) g ~ 

r fuxffl / v^(/(i-^)) \ , ( uxV -( H ')) 7 (i + « 2 ) / ( t . 

7 '" (nx MTmT (I - ^) — v M 5o °( h X H ; 



r " = " fidfe) + 2 (( 2 ^ r '+ r o-5oo5°T°) x u) 



where we used vectorial notations (h°) Q = h° a , (h l ) Q = /i^, as well as the vectors (H 1 )^ = Hp — gonh° (h 
the following various vectors (ff) = ff , (ff)^ = ff , (r?)' 3 = rf , and (rj])' 3 = rjf '. The matrix J is given by 
and the effective mass by fa = mgoo^o ( i-u* ) as we ^ as 9* = 9oo9° l - The expression Eq.([T4| for 



J = 

the energy operator H will be the one to diagonalize in the next section. 



SEMI-CLASSICAL ENERGY 



The semiclassical diagonalization the Hamiltonian Eq. (fl4l) is expected to lead to an effective Hamiltonian with 
gauge fields resulting from the back reaction of the spin degree of freedom (fast) on the translational momentum which 
can be treated semiclassically for slowly varying enough inhomogeneities. Indeed, the emergence of gauge fields is a 
general feature of systems providing fast and slow degrees of freedom. The purpose of ref [9|, was to investigate the 
origin of quantum gauge fields and forces by considering the diagonalization of an arbitrary matrix valued quantum 
Hamiltonian. To be precise, by diagonalization we mean the derivation of an effective in-band Hamiltonian made of 
block-diagonal energy subspaces. This approach, based on a new differential calculus on a non-commutative space 
where h plays the role of running parameter, leads to an in-band energy operator that can be obtained systematically 
up to arbitrary order in H. Particularly important for our purpose, it has been possible, for an arbitrary Hamiltonian 
H(R, P) with the canonical coordinates and momentum = iH5?, to obtain the corresponding diagonal 

representation e (r, p) to order h 2 , in terms of non-canonical coordinates and momentum (r, p) defined later and 
commutators between gauge fields. The method is quite involved, so that in the present paper we restrict ourself to 
the semiclassical approximation (order h). 

The mathematical difficulty in performing the diagonalization of H comes from the intricate entanglement of 
noncommuting operators due to the canonical relation = ihSf. In Q starting with a very general but time 

independent H (R, P) and by considering h as a running parameter, we related the in-band Hamiltonian VHV + = 
e (X) and the unitary transforming matrix V (X) (where X = (R, P)) to their classical expressions through integro- 
differential operators, i.e. £ (X) = O (e(Xo)) and V (X) = AT(y(Xo)), where in the matrices e(Xo) and V(Xo), the 
dynamical operators X are replaced by classical commuting variables X = (Ro, Po)- 

The only requirement of the method is therefore the knowledge of V(Xo) which gives the diagonal form e(X ). 
Generally, these equations do not allow to find directly e (X), V (X), however, they allow us to produce the solution 
recursively in a series expansion in h. With this assumption that both e and V can be expanded in power series of H, 
we determined in Q, the explicit n-th in band energy to order h for an arbitrary given Hamiltonian 

e n (r, p) = e ,n (r, p) + jV n { [e (r, p) , A Rl ] A p ' - [e (r, p) , A Pl ] A R ' } + 0(h 2 ) (15) 

where Ar = ihW-pV + and A-p = — ihW-RV + with V the diagonalizing matrix. The operator V n has the meaning 
of the projection on the n-th energy subspace. The new non-canonical dynamical operators r and p depend on gauge 
fields An = V n {Aa^) and Ap—VniAp) similarly to electromagnetism, as we have r = ihd p + An and p = P + Ap. 
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These gauge invariant quantities not only are emerging naturally but are also necessary to have a gauge invariant 
energy Eq. (fl"5|) . The operator eo (r,p) is the diagonal energy obtained at zero order (h°) in which the classical 
variable R and Po are replaced by new non-commuting operators r, p. Therfore the first step consists in finding the 
matrix V(X.q) which diagonalizes the "classical" Hamiltonian H(K ). 



Zero Order Diagonalization. 

For practical purpose we introduce the three dimensional effective metric G 1 * — H l a Hp8 a ^ as well as the gravity 

coupled momentum P a =H l a (Pi + |e e/ 3 7 f f'V 7 ). 

As shown in Q , the classical block-diagonalization of the Hamiltonian Eq. (J14I) , but without the last term (cor- 
responding to the static case), can be performed by the following unitary FW-like matrix (denoted Fq for Foldy- 
Wouthuysen) 



F (P) = D ( E + m + c-(3 [a.P + P .or) + iV ) /^/2E Q (E + m) 



* = /("^J N - and D = I + fr gifl^ffi* 

Indeed one can easily check that 



F HF + = pjPiG»Pj + he^T^^G^P + m 



! ) — 



(r° + r e ). mE + i- ' 1 

' ( 300 (h° x h*) u x H l 



,r \" " / ■ mS+ ^ ' + -g*P i + -P i9 * 

4^0 V 7 1 [Eq + m) / 2 2 



(1 -u 2 ) 



(Viu) J F + (16) 



with fi 7 = EgfjjTf^ . All contributions are block-diagonal except the last term which was not present before in the 
case of a static metrics 0) ■ The proof of this block diagonalization relies on the simple fact that for classical variables 
Xq, the matrices H l a and f" are independent of both the momentum and position, f3 and a.P anticommute and in 
the Taylor expansion of Eq all terms commute with j3 and a.P + P+.a. 

As said before, the last term in Eq.tfTfJf is non diagonal and must treated specifically. Actuallt, one can apply a 

second unitary transformation that will cancel the non diagonal contributions of Fq( ( g00 ( h J*^)"* 1 * - . (V,u) J)Fq + 
without affecting the rest of the diagonalized Hamiltonian to the first order in h. The explicit form of this transfor- 
mation is : 

_ 1 _ v _ / fan (h° x h') u x H') JF<>+ \ „ 



(1-u 2 ) v ' J 2^/P^Pj+m 2 

V- being the projection outside the diagonal. One can check that Fq — 1 is antihermitian so that Fq is unitary to 
the first order. As a consequence, the composition of the two unitary transformations yields the following diagonal 
energy operator : 



e (R,P,i) = F^FqHqFq+F^ = pJPiGvPj + he aM Tf Yn&iPi + m 2 



( r + r ') • + wM) ) + l s ' p < + \™ 
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The last term is explicitely given by: 

' ( 5 oo (h° x h*) - u x H* 



V+ F 



(1-u 2 ) 



(Viu) J F H 



DV+ [E + m + c/3 a.P 



(E + V(r)m) J + c/3 (e.p)) ^ M ^ + 



= ^ (ggo (h°xhO-uxff).(V iU ) - s _^ ((/- 1 )gooj h ° xh y-( V ^) p, s 



2^o (1 - u 2 ) 
Ultimately, the diagonalization process yields: 
e (R,P,t) = F{ t F () H () F a + F{+ 



2fE Q (1 - u 2 ) 



(f + r e ) 



TO.E 



a 1 Pi 



(E + m) 

where for the moment R and P arc treated as classical commuting quantities 



2 ,P,+ 2 P,S+ ^ 2/£„(l-u 2 ) RS 



(18) 



First order in h diagonalization 

From expression Eq. (|18|) we can deduce the diagonal energy operator for, let say, the particule subspace £q,+- The 
semiclassical energy is given by Eq. (fT5|) . where £o,+ corresponds to the positive energy subspace Eq. (fl8|) in which the 
classical variables R, P are replaced by the quantum covariant ones r = ihd p + HAr and p = P + hAp. The explicit 
computation for the Berry connections Ap = V+(Ar) and Ap = V+(Ap) with An, = V+(ih {FqF ) X?p(FqF ) + ) an d 
Ap = V + (ih(F/ ) F )V R (F^F ) + ) gives the components 



2E (E + to) 



ik + o (h) 



(19) 



-4, 



2E (E + to) 



o(h) 



(20) 



where E is the same as Eq above, but now R is an operator and is the inverse matrix of Hp. 

We also denote, for the rest of the paper, p to be the same expression as P in which R and P have been replaced 
by r and p, namely: 

p a =H i a { Pi + \E eM ff<T') 

To complete the diagonalization we need to evaluate the quantity 

M+ = \ v + { h (x) ' A * l \ ^ h (x) ' A °] A * 1 } + ° {h2) 

which being on all point similar to the one given in Q or Q is simply stated : 

M+ = i (Qe - (A R xp)) .B - ivm(r). (p x E) 

where the " magnetotorsion field" B is defined through a three dimensional effective torsion tensor B 1 
— ^PsT Sa Pe a p 1 where the effective torsion T Sa/3 is defined as 



Ta ps = jgs( H la d,H k0 - H lf3 d,H ka ) + H la ff 5 - H w tf 
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The physical origin of the term M has been discussed in Q and Q for the static case. Note that for the static case 
( g ° l = 0) we retrieve the true torsion T a|55 = H s k (H la diH k(s - H lf} diH ka ) + H la T^ S ~ H w Tf s . 

Considering similarly the anti-particle subspace, it turns out that the coordinate operators are identical with the 
particule ones and that M_ = — M + therefore the full energy operator with both particles and anti particles can be 
cast in the form 

1 * 1 - » ft* , r e\ , ( S -H^)^P^ ^((/-l) 5oo (hOxhO).(V 4 u)^ 

(21) 



e (P, r) = 0+-^+-^+- ^ + r*j . frill + ^ +/W ^ 2fE(l-u^ ^ m 



where M = M + and 



4£ ' \ {E + rn)JJ y AE \ (E + fh) 

with the vector f , defined in terms of its components as (pi) = (r). Equation (f2"Tj) is the main result of 

this paper. 

The coordinates operators in Eq. (|21l) satisfy a non commutative algebra as 

[r i ,r j ] = ih 2 e% (22) 

\p hPj ] = m 2 e% (23) 

[Pi,rj] = -ih 9ij +m 2 &% (24) 



where S^ 8 = d ai Ap j — d/3 j A oli + [Aq^.A^J the so-called Berry curvature. An explicit computation gives : 
QPP = _Ji I mS, 1 V y 



i / .„ . ^ ^; p 7 \ ~ ~ „ fl „ i 



« 2P 3 1 (J5 + m) 

efj = ^ (mE 7 + ) ^V n P a ^, - ^V ri m (E x P) (25) 



There are also other Berry curvature mixing coordinates and spin 



< = I*. ^ = -^ "^% + + P ^ J ^ Vr ^ (26) 



Interestingly Eq. (|2"Tj) can be rewritten as 



£ = (pi - & (pi - |f + m»(r) + i 9 V,: + ^,9 - J (f + T e ) .0" 

+ Phc^ TfE(l-J) P-S+^M (27) 



where 0^ r = ^ ^?™£ 7 + (e+L) " ) * s "rescaled" Berry curvature. This formula clearly shows that the spin 
connection couples only to the Berry curvature. 

We note in Eq. (j2"7) the presence of the term — | + .0 rr + \g % p% + \pig l not proportional to j3 and 
independant of the particle charge, that is discrimnating between particles and anti particles. These terms give 
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different energy levels for the Dirac particles and antiparticles as e+ 7^ — e_. The coupling 
proportional to the spin will survive in the case of the non-diagonal static gravitational field but cancels for a diagonal 
metrics as studied in [2j. On the other hand the term \g l pi + \pig l vanishes for a static metrics since in this case 
g 1 = d g lJ = so that g 1 = 0. Therefore the symmetry between particle and antiparticle is restablished only for 
static diagonal metrics. 



THE STATIC GRAVITATIONAL FIELD 



This case is caracterized by the following time independent metric: g^ = gij (R), 500 = .9oo (R)j 9i0 — 0. In that 
case expressions simplify greatly. Actually, the transformation matrix U is U = y/—gh,Q and = v/ yf—gh^ = f. 
The effective quantities reduce to: 



(h 4 ^ = Hp = Hp = goohohp — y/goohp 
Y-qP _ _ pe/3 

i i i 

r° = f ° = \ [e^Tf + e^Hprr} = -~<?oo (fcg) 2 e^Wph^Tf 

r e = 

G ij = //,,,//" 
m = ?ri(?ooft.Q as w = 

where r*' stands for the Christoffel symbol. In this case Eq. (|2ip reduces to 



e * ( ft + A Fi (r) . (rnE + |||)) & (r) ( w + ^r, (,) . (™E + |||)) + r& + 

+ Af° f^ E , + ( S -P)j^ (28) 
4E 7 V (E + m)J y ' 

with p a =y/gooh l a {pi + j £ e /3 7 r^cr 7 ) and the vector I\ defined in terms of its components as (1^) = e Q ^ 7 r"^ (r) . 
Then even for this case, particles and antiparticles avec a different in-band energy operator because of term 
jgf* (mT, 1 + ^fiJyj breaking the symmetry e + 7^ — e_. Although this term was already in previous studies [HQ 

this essential point has not been pointed out. For a diagonal metric To = 0, and the symmetry particles/anti-particles 
is recovered. 



ULTRARELATIVISTIC LIMIT 



It is interesting to look at the ultrarelativistic limit mc 2 —t 0. One readily obtain 

^ ((/-l) g oo(h"xhQ).(V t u) Xg 00 Bp A ^ft+r 7 ) 1 1 

6 " ^ + ^ 2/(1-11*) A + ^ WT" 4 I + 2 3 ^ + 2 Pi9 (29) 



with e = cy ^pj + f P ) G 4 - 7 + ^ P ^ and A = Hp.'S/p a biased hclicity, that is not projected on the 

momentum p but rather on p. This fact is not astonishing. Actually, as we saw in the diagonalization process, the 
particle is submitted to the action of an effective gravitationnal field, which differs slightly from the initial field. This 
effective metric is responsible for considering the momentum p rather than as a dynamical variable P. Nicely this 
energy can be expressed in terms of the helicity which is the relevant variable for massless particles and not in term 
of S. As shown in [(|, Eq. (|29|) is also valid for photon with the one-half spin matrix £ replaced with spin one matrix 
S. 

Here also we see that photons and anti-photons do not have the same energy spectrum. The symmetry is again 
only restored for a static diagonal metric. 
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THE TIME DEPENDENT SYMMETRIC GRAVITATIONAL FIELD 

A typical example of such a metric is the Schwarzschild space-time in isotropic coordinates. This case, studied 
in a different manner in 0and for a time independent metric, received a full independent treatment within our 
formalism in We now present the equivalent results for a time dependent metric, completed with the spin 

matrix dynamics. For a symmetric metric, the semiclassical Hamiltonian has the following form [3[ 

H = ]- (a.PF(R, t) + F(R, t)a.P) + f3mV(R, t) (30) 

corresponding to the metric gij — Sij ( ^-pjf£^ ) 7 9m — and 500 — V 2 (R,t). We will define <f) = j,. In that context 
the relevant quantities for the diagonalization appear to be : 



and 



yeP _ pe/3 _ ygp 
i i i 

0, J°4>M 



r 



4>v 



^_ f0 _h F (R) p F(R) f fV(R)\\ h i h »-() 

r e = 

Qij = y2 g ij 

Similar computations to the ones performed in the previous section lead to the following expressions for the dynamical 
variables and the diagonalized Hamiltonian 

r = R-/i ±- LJ p = p (31 

2E(E + mV(R)) v v ' 



and the diagonal energy becomes: 



F 3 (r,i) 



e = (3JF 2 (r,t)P 2 + P 2 F 2 (r,t) + mV 2 (r,t) -^-mhpV 4>(r,t). (P x S) (32) 

The Berry curvatures are given by: 

% = _^M mV ^(r) (S x P) . (34) 

eg - " (35) 

and Q r z f being unchanged, Q$f = 0, and e = y/F 2 (r,t)P 2 + P 2 F 2 {v 1 t) + mV 2 (r,t). 

Note here that the magnetotorsion field B = 0. From Appendix B, we can also get the non Hcrmitian contributions 
the time evolution operator which in this case reads jftjj In / = ^ lny'— gh®. 

One can check, after developing r as a function of R and the Berry phase, that our Hamiltonian coincides, in 
the weak field approximation,with the one given in Q when considering the semiclassical limit (order h). This also 
confirms the validity of the Foldy Wouthuysen approach asserted in Q . 
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ONLY TIME DEPENDENT METRIC 



The metric tensor and the vierbein only depend on time. Therefore we have 



(H 1 )^ = H} 3 = g 00 h I h) } ^- 



(*), 



7 / I TJ% 

tl a — tin 



2g ao (h ^ - ftO/ij) u 5 



(1-u 2 ) 



r° 7 



ff = rf - fe^ 7 {H-% gm h^ K ^ 



r; = o 

so that the energy becomes 

fi / (s.H»p,)H« K \ ! . ! 

e (p, r,t) = ^+^r°. mS + V + -g*(t) Pj + - Pi g\t)+phM (36) 

with 

. - cfc + Ar. ,,«, . (fcs + pi)) G ,. m (» + ^r. ,,«) . (*= + ff| 

which show that particles and antiparticles have a different energy spectrum in this graviational field. 



CONCLUSION 



The semiclassical limit for Dirac particles interacting with a fully general gravitational field was investigated through 
a first order in h diagonalization of the Dirac Hamiltonian. This work extends previous ones where only static metrics 
were considered. The time dependence of the metrics leads to new contributions of the in-band energy operator. In 
particular we found a coupling term between the linear momentum and the spin, and terms which in general will 
break the particle - antiparticle symmetry. 

As already found by other authors, the time dependence leads also to special features like the non-unitarity of the 
evolution operator, whose origin can be tracked back to the notion of scalar product in the Hilbert space of wave 
functions for a time dependent metric. This non-unitarity is unavoidable but we could nevertheless diagonalize the 
full evolution operator, even though our main focus was to obtain the block-diagonal form of the energy, this one 
turning out to be Hcrmitian. In addition, to the very general semiclassical diagonal energy operator, we provided 
several physically relevant examples. 



[1] L. L. Foldy and S. A. Woutfmysen, Phys. Rev. 78 (1950) 29. 

[2] Y. N. Obukhov, Phys. Rev. Lett 86 (2001) 192; Fortschr. Phys. 50 (2002) 711; Phys. Rev. Lett 89 (2002) 068903. 
[3] A. J. Silenko and O. V. Teryaev, Phys. Rev. D 71 (2005) 064016. 

[4] J. D. Bjorken and S. D. Drell, Relativistic Quantum Fields, New York: McGraw- Hill, (1965). 
[5] M. Leclerc, Class. Quant. Grav. 23 (2006) 4013. 



12 



[6] P. Gosselin, A. Berard, H. Mohrbach, Phys. Rev. D, 75 (2007) 084035. 
[7] P. Gosselin, A. Berard, H. Mohrbach, Phys. Lett. A, 368 (2007) 356. 

[8] P. Gosselin and H. Mohrbach, Eur. Phys. J. C 64, (2009) 495, P. Gosselin, J. Hanssen and H. Mohrbach, Phys. Rev. D 77 

(2008) 085008, P. Gosselin, A. Berard and H. Mohrbach, Eur. Phys. J. B 58 (2007) 137. 
[9] P. Gosselin and H. Mohrbach, J. Phys. A: Math. Theor. 43 (2010) 354025. 



APPENDIX. 



Derivation of Eq. ()14[) . 



Let us start with the following development for H: 



H = goohPrf-fP* + \e ePl Tf^ + i|r°V + goohPrfm + \d t In (-gg oa ) 

= g 00 g Oi (Pi + he^^v + th^a-y) 

+ z.g 00 (h° s h^ K ) S K (F 4 + he et)l ^Y? + ih^cP) 

pg/3 p07 

+ g m hl^h l {P % + he eM ^-TP + ih-^of) 

■pg/3 p07 

- gooofih^Pi + hE e p^JP + ih-±-cP) 

+ \e eM Tf^ + i^faf + ftohtfm + 5ft In (-99°°) 
and remark that we can rewrite the four first terms in the following form : 



3oo 



/ h a h i \ v el3 r° 7 



pg/3 p07 



9oo9 0i (Pi + he efh -L-ir + ih-±-aP) 
- ig o {hPsh^e^) £ K (P 4 + Ke e01 ^-YP) 

p0 7 

■ 9oohW U - 



x\p + ne eP ^ - ( - (_| m ^ 2^ + in^oP 



reP p°7 

9oog 0i (Pi + fe^-^s 7 + 



T el3 



Woo 



(ftXe^)£ K (P < + fe e/37 4-5T) 



/ pg^ p07 \ pe/3 p07 

= a. Hi Pi + he etil ^-YP + ih^-cP + g 00 g m (P + he^^-YP + ih^-cP) 



4 4/ v ' yy ' 4 4 
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where the effective dreibein and spin connection : 

of i h l h o\ 
Hp = gooK I hp j^p- I 

r o 7 

f f = Tf - (H-^ g 00 (h° s h^ K ) -f 

have the form claimed in the text. 

Now, we compute H + H + , which is the first contribution in Eq.© (up to the factors jjrz^^w) that will be skipped 
constantly in this section for the sake of readability, and reintroduced ultimately) . 



1 / \ 1 / t efs \ 1 / V ef} \ 

+ \goog 0l P t + l^goog 01 

+ goo gM(He eM -^W) - hV Rt (goo {h° s h^)) S K - g 00 {h Q 5 h^ K ) ^"(te^-i-jE,, 
4 4 



1 ( „ , tf\ 1 /„ , f eP 



2 o^IPp I P< + he^-^-JPj + - I Pi : + Ke e01 ^Yn\ a fi 'Hp + g 00 h° j3m 
+ (f - RVj* (.900 (ftgft*e%))) .E + lgoo9 0i Pi + \Pi9m9 0i 



with 



) = ^e^7 r o + hg 00 g m £ e p 1 -^ 



pe/3 po^ 



Introduce, as in the text (r )^ = T 07 = e e pyTg^, (r 1 )^ = goog^hSgp^Tf 13 , and use the following expression for the 
spin connection 

vf=h%h$ (ftfir-vft?) 

to show that %oo {hghts VK ) £ k ^\£q@u^4-) = 0. As a consequence, we are left with 

f = ^(r + fl0off w r i ) + RH i x^ 

as announced. 

To compute the other contributions to U^HU^? in Eq. (|SJ), we need to find the expressions for 
—upa^ ^ H+ ^ - Upa^, ' M<3 " - , and _ ) v p-^J^( m ' 9 " ) ^ — upa^. We start with the computa- 

tion of upa 13 (h + upa^ by decomposing (^H + H + ^j as: 



H + H+) [H 1 + H 



T'+ 



2 



r° - ftV ft ( 500 (h° xh 1 )) EH g 00 h° j3 



lgoo9 0i P + \Pi9ao9 m 



where 



(H' + P'+) 



1 / f e/3 \ 1 / t el3 \ 

-a^H} [ P + fe^ 7 ^£ 7 ) + - ( Pi + fe^ 7 ^£ 7 ) ofiH} 



. / fe/3 \ . ~ 

For the sake of the computations, we will denote P a = H l a [Pi + fe e/ 3 7 ^-S 7 J = H^Pi.so that 

(H> + H'+) i i 

'-una! 5 = -upo p ofiH % p P t upofi + -upoP P l of upoP 



iu 7 a 7 a f3 Hj j P t u 1 a' y + ^u^a 1 P l a^ 3 u^cf Hp 



1 

2 
1 

2 
1 

2^7 



1 



-^/goou^a' 1 a p u^a" 1 Hj 3 P l + -P % HJ s u 1 a' 1 'or u 1 a 





1 




P 4 ,u 7 a 7 


+ 2 


u 7 a 7 ,P Z 



7 oA 



1 



— a^cP oP u-yCp HpP 1 + —P l HpUyCt ' n (/- n 



7/35 Qs«flj [p\u 7 cr 



P\u 7 a 7 



1 

2 

= ^a' 3 u 2 HlP l + ip'iPguV 
+ iu^^u^H^P 1 - HpP'u^iu^^H 5 



+ -iUjE 



7/35 



E 6 ^ 



P\ u 7 a 7 



P 4 ,u 7 a: 7 



+ 



(u 2 £^ - (uS) u") P* P 4 - P*,P J (u 2 £^ - (uS) u") 



Now, given that: 



7/35 



P l ,w 7 a 7 



-itt 7 e 



7/35 



S 5 ^ 



P , u 7 a 7 



P\w 7 a 7 



P J , u 7 a 7 



-ViU 7 a 7 



\Q0 



-hu^Hl (-2 JV^ + f f u c ) 



one thus has, 



[Pi' + 



\oPm 2 P^ + \ppu 2 af 3 - ^g 00 h° hu^ K5 Hi (-2 JV t u s ) 



since, by an argument already used, u 7 e 7Kl5 P*f f 5 



■ . v ih ' x h^ff u e = 0. 
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To complete the computation of upa^- — g — -up<y? , we need to calculate the following contribution: 
u^of (g 00 h° Pm + i 5o o5° 1 ^ + ^W* + (f ° - hV Ri (goo (h° x h*)) 7 ) u 7 a 7 
= -. 9 oo^/3u 2 m + ^g og 0i u 2 Pi + ^u 2 P igoo g m + u 2 (f° - ftV^ ( 5o o (h° x h 4 ) 7 )) E 7 
- 2u 2 (f° - hV Ri [goo (h° x h^)) S 7 + u.S (f° - SV* ( 5o0 (h° x h') 7 )) u 7 
= -goohlPu 2 m + -ffoo^u 2 ^ + iu 2 p 5oo5 w + u.S (f° - (.g 00 (h° x h<) 7 )) « 7 
-u 2 (fO-ftV^ (.goo (h°xh J ) 7 ))s 7 

Now, we turn our attention toward the second contribution in Eq. ©, namely [^' Uf3Cl ' M<3 ° 1 . To do so, and 



since 



H, upa^ 



1 r 
2 



we only need to concentrate on the anti hermitian part of H. Given that, 



we can write that 

ff-iP 1 



and 



-1O7 



-1Q7 



H = a? Hi P + fe e/3 * + iK^-cC 1 + goog 0l ih^-c 



igoo (h° s hi,e s ^) £ K (P + he eM ^-V) + i-r°V 



^0/3 



-,03 



/3e 



P, u^a^ 



H+.upa 13 



F o/3 r o,3 



-ft I V 4 (4) + .9005°'^ + -f - - H^j 

- .900ft (h^e 5 ^) JV iUk 
Ultimately, we need the third contribution to Eq. (jHJ: 



l 2 



1 — 
2 



UpOi' 



! , (v p p) .V R 



-KE. (u x ET) 



[/(l-t^)] 2 

vg. (/(i-" 2 )A 

[/(1-u 2 )] 2 / 



1 * rT 1 



We can now gather all these terms to obtain the expression of U 2 HU 
Define, as in text, the vectors Tf, T° by: 

,/3 



(ff) =ff,(r°)" = r^(rg)^ = C 
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and H%: 



(h° s hi - h$hi) us 



TTl TT'L 



Ultimately, reintroducing the factors nJ"^) when needed in Eq. §E§ , U 2 HU 2 can be written in a compact form 



n (l + U 2 ) 1 n . 1 

+ gooK ^ _ u2 j (3m + -goog^Pi + ^Pigoog 



(l» + r-) D - ,fa.(-°xhW 8 uxH-) (Vfu) ; 



where we have defined T° by 



.„ (i + u 2 )«-u^. n // ... rR n r° 

->0 V / 7 7-pO , t [ [ ffipJ 1 M) _ „0i z 

7 ~ (1 _ u 2) i v + n y[ H 5 T i +Y~ 9009 T 



7 



r e 

7 



fuxH , / V^(/(1-U^)) \ (uxV < (HQ) 7 (l + u 2 ^^ , f h o xh ^ 
(u x H )^ ^ /2(1 _ u2) 3 ) + h (1 _ u2) - ft (I" ^ V ^ (^00 (h x h )J 



Non Hermitian contributions of the time evolution operator 

We compute here the contributions to the diagonalization due to the non- hermitian term —ihU? (t) ^ (£) 
First we obain 



t«7» (i) -c/ (<) - ih /(1 _ u2) + ^ (1 _ u2) + * fi (T-Ga 



TP ^t_J — 



(l-u») 



(37) 



and then, the contributions of the terms in Eq. (|37p to the diagonalization are computed by applying the transfor- 
mation matrix and then projecting on the diagonal blocks. We are left with: 



U 



~\ + 



u[p) =ih 9t L ;. y ~ — jr^ + zft- 

/ (1 - u 2 



(1-u 2 ) 



ihc/3 — m 



(l-u 2 )£ 



He 2 (Px (Px (ux|u 



